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Scaling  Phenomena  ? 


•  Detection:  scaling  ?  what  does  it  mean  ?  non  stationarity  ? 

•  Identification:  relevant  Stochastic  models  ? 

•  Estimation:  relevant  parameter  estimation  ? 


•  Side  Issues: 

Robustness  ?  Computational  Cost  ?  Real  Time  ?  On  Line  ? 
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Outline 


I.  Intuitions,  models,  tools 

1.1  Intuitions,  Definition, Applications 

1.2  Stochastic  Models:  Self-Similarity  vs  MultiFractal 

1.3  Multiresolution  Tools,  Aggregation,  Increments 

1.4  Wavelets,  Continuous,  Discrete 

II.  Second  Order  Analysis,  Self  similarity  and  Long  Memory 

11.1  Random  Waks,  Self  Similarity,  Long  Memory, 

11.2  2nd  Order  Wavelet  Statistical  Analysis, 

11.3  Estimation,  Estimation  Performance, 

11.4  Robustness  against  Non  Stationarities, 

III.  Higher  Order  Analysis,  Multifractal  Processes 

111.1  Multiplicative  cascades,  Multifractal  processes, 

111. 2  Higher  Order  Wavelet  Statistical  Analysis, 

111. 3  Finiteness  of  Moments, 

111. 4  Estimation,  Estimation  Performance, 

111. 5  Negative  Orders, 

1 1 1.6  Beyond  Power  Laws. 
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Irregularities,  Variabilities 


Scaling  or  Non  Stationarities? 


[3] 


nb  connexions 


Scaling  ? 


Trafic  (WAN)  Internet 


Trafic  (WAN)  Internet 
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Scaling  ? 
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Log  (DSP) - Nombre  Octets 


Scaling  ? 


Trafic  (LAN)  Ethernet - Densite  Spectrale  de  Puissance 
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nb  connexions  nb  connexions 


Scaling  ! 


•  Definition  : 

Non  Property:  No  Characteristic  Scale. 

Non  Gaussian,  Non  Stationary,  Non  Linear 

•  Evidence: 

The  whole  resembles  to  its  part,  the  part  resembles  to  the  whole. 


•  ANALYSIS:  Rather  than  for  a  characteristic  scale, 
look  for  a  relation,  a  mecanism,  a  cascade  between  scales. 
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Scaling  :  Operational  Definitions 


•  Multiresolution  quantity: 

Tx(a,  t)  (e.g.,  Wavelet  Coef.). 


•  Power  Laws: 


i 

n 


E|  Tx(a,t)\q 
J2nk=1\Tx(a,tk)\ q 


—  Cq 
=  C 


Q 


a 

a 


C  (<?)i 
C  (q) 


-  FOR  A  RANGE  OF  SCALES  a, 

-  FOR  A  RANGE  OF  ORDERS  q 

-  Scaling  Exponents  ((q). 


•  Beyond  Power  Laws  :  Warped  Inf.  Div.  Cascades 

E|Tx(a,  t)\q  =  Cq\a\^q^  =  Cq  exp(£(g)  In  a) 
E|Tx(a,t)|9  =  =  Cqexp(C(q)n(a)) 

— >  Visit  Pierre  Chainais’s  Poster 
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Ubiquity  ? 


Trafic  (LAN)  Ethernet  Trafic  (WAN)  Internet 


temps(s)  temps(s) 
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Ubiquity  ! 


-  Hydrodynamic  Turbulence, 

-  Physiology,  Biological  Rythms  (Heart  beat,  walk), 

-  Geophysics  (Faults  Repartition,  Earthquakes), 

-  Hydrology  (Water  Levels), 

-  Statistical  Physics  (Long  Range  Interactions), 

-  Thermal  Noises  (semi-conductors), 

-  Information  Flux  on  Networks,  Computer  Network  Traffic, 

-  Population  Repartition  (local:  cities,  global:  continent), 

-  Financial  Markets  (Daily  returns,  Volatily,  Currencies  Exchange  Rates), 
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Analysing  Tool  1  :  Aggregation 


Compare  Data  against  a  Box,  then  vary  a 

Tx(a,t)=^Jtt+aToX(u)du 

Average 


1/T 

- . 

T 

-10  12 

temps 

Works  only  for  Positive  Time  Series,  Density 


ini 


Analysing  Tool  2  :  Increments 


Compare  Data  against  a  difference  of  Delta  Functions,  then  vary  a 

Tx(a,t)  —  X(t  +  aro)  -  X(t) 

Difference 


-10  12 
temps 

Increments  of  higher  orders  or  generalised  ^-variations 

-  Order  2  :  Tx(a,t)  =  -X(t  +  2 ar0)  +  2 X(t  +  ar0 )  —  X(t), 

-  Order  N  :  Tx(a,t )  =  S^Lo(_1)Pap^(t  +Paro), 
where  E^=o(_1)PapPfc  =  k  =  0,  •  •  •  > N  -  1. 
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Analysing  Tool:  Multiresolution  Analysis 


•  Multiresolution  Quantities: 

X(t)  -  TX(a,t)  =  (fa,t\X),  faM  =  i/o(V) 


Aggregation 

fo(u)  =  (/30) 

=  ^kl!+aT°x(u)du 

Box,  Average 

- 1 - ! - ! - 


1/T 


T 


temps 


Increments 

fo(u)  =  (io) 

=  X(t  +  aro)  —  X(t) 

Difference 


-10  12 
temps 


? 

? 

? 

? 


? 
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Analysing  Tool:  Multiresolution  Analysis 


•  Multiresolution  Quantities: 

X(t)^Tx(a,t)  =  {fa>t \X), 

•  Choices  for  Mother  Functions: 

Aggregation 

fo{u)  =  ((3q)*n 


=  It+aT°  X('u')du 

Box,  Average 


1/T 

- 

- 

T 

- 

temps 

N 


fo, 

Increments 

fo(u)  =  (. Io)*N 
—  X(t  +  cltq)  —  X(t) 

Difference 


-10  12 
temps 

N 


/o(^) 

Wavelets 

fo(u)  =  1pO,N 

=  fX(u)  ^o(^), 

Average,  Difference 


-10  12 
Temps 

N 
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Wavelets  and  Scaling:  Key  Ingredients 


Dilation  Operator,  -^o(-£-) 


Temps 


•  Number  of  Vanishing  Moments, 

N  >  1,  f  =  0,  k  =  0, 1, . . . ,  N  —  1. 
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Wavelet  Transforms 


•  Mother-Wavelet  and  "Basis”:  /  ip0(u)du  =  o,  ipa,t(u)  =  ra^o(^) 

•  Wavelet  Coefficients:  continuous  wt  Tx(a,  t)  =  (X,^a>t) 


SCALE 


TIME 


Modulus  Maxima  WT 

Skeleton  :  Maxima  Lines 


and  Discrete  WT 

dx(j ,  fe)  =  Tx(a  =  2-y ,  t  =  2Jfc) 


WTMM 


SCALE 


TIME 
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Outline 
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Mod.  Tool  1 :  Rand.  Walks  and  Self  Similarity 


Random  Walk:  x(t  +  r)  =  x(t)+  sTx(t) 

Steps  or  Increments 

Statistical  properties  of  the  steps: 

-  A1 :  Stationary, 

-  A2:  Independent, 

-  A3:  Gaussian, 

=>  Ordinary  Random  Walk,  Ordinary  Brownian  Motion, 
=>  E X(t)2  =  2D\t\,  Einstein  relation, 

=*  E|X(i)|«  =  2D\t\^2,  q  >  -1. 

Anomalies: 

=*  E X{t)2  =  2D\ty, 

=>  E X(t)2  =  oo. 

Self  Similar  Random  Walks: 

-  B1 :  Stationary, 

-  B2:  Self  Similarity 


[18] 


Modelling  Tool  1 :  Self-Similarity 


•  Definition:  STX(t)  f=  cH5T/cX(t/c ),  Vc  >  0,  Dilation  Factor, 

l  >  H  >  0  :  Self-Similarity  Exponent 


•  Interpretations: 


Covariance  under  Dilation  (Change  of  Sc^le}, 
The  Whole  and  the  SubPart  (Statistically) 

No  Characteristic  Scale  of  Time. 


v 


B 


a>in/-\iip 

t,  i  ,,  lOl 

rnj  m 


v 


LE, 


Implications: 

-  Non  Stationarity  process  with  stationary  increments 

-  IE| X(t  +  a r0)  -  X(t)  lq  =  Cq\a\qH , 

-  Va  >  0,Vc  >  0,Vg  /\E\X(t)\q  <  oo, 


-  a  Single  Scaling  Exponent  H. 

-  Additive  Structure, 

-  (Correlated)  Random  Walk,  Long  Memory,  Long  Range  Correlations. 
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Mod.  Tool  1  (bis):  Long  Range  Dependence 


•  Definitions  : 


-  Let  X  be  a  2nd  stationary  process  with, 

-  Covariance  :  cx(j )  =  IE x(t)X(t  +  r) 

-  Spectrum  :  rx(v) 


CX(t) 

=  cT|r|  0, 

0  <  (3  <  1, 

r  — >  +oo 

=  c/M_a> 

0  <  a  <  1, 

M  -►  o 

WITH  a  =  1  —  (3  AND  cf  = 

2(27r)  sin((l  — 

7)tt/2)ct  . 

•  Consequences  : 

■  J2a°°  cx(T )dr  =  +  oo,  A  >  0, 

-  No  Characteristic  Scale, 

-  Aggregation  :  Tx(a ,  t)  =  ^  Jtt+aT°  X(u)du, 

=>  VAR  Tx(a,  t)  ~  Caa_1,  a  — ■>  +oo, 

-  Increments  of  Self.-Sim.  Proc.  (with  H  >  1/2) 
are  Long  Range  Dep.  (with  a  =  2 H  -  l). 
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Wavelets  and  Self-Similar  Processes  with 
Stationary  Increments  -  Summary 

(Flandrin  et  al.,  Tewfik  and  Kim) 

•  PI :  {dx(j,  k),  k  e  2}  Stationary  Sequences  for  each  Scale  2j  . 

N  >  1 

•  P2:  Self-Similarity  :  Dilation 

{X(t)}  =  {cHX(t/c)}  =*  {dx(0,k)}  =  {2 -jHdx(j,k)} 

-  P3  :  Marginal  Dist.  P:i(d)  =  fto  =  ffrj  . 

•  P4  :  {dx(j,  k)}  Short  Range  Dependent  if  n  >  H  + 1/2. 

|2Jfc  —  2J  k'\  — >  -poo,  \Cov dx{j,  k)dx{j' ,l)\  <  D\2^k  —  2J  k'\2<yH~N\ 

N  >  1  and  Dilation 
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Daubechies2  Haar 


Wavelets  and  Long  Range  Dependence 


(Flandrin) 


V 

Uwv  p,vy  Y'v 

Nyv  A/v^yirv  wyyy1 

v‘Yvy  w\jy  vyu 

uv^wl/Wy  Wv 
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Wavelets  and  Self-Similar  Processes  with 
Stationary  Increments  -  Summary 


•  PI :  {dx{j,  k),k  g  2}  Stationary  Sequences  for  each  Scale  2j . 

N  >  1 

•  P2:  Self-Similarity  :  Dilation 

{X(t)}  =  {cHX{t/c)}  =*  {dx(0,k)}  =  {2 -iHdx(J,k)} 

-  P3  :  Marginal  Dist.  P:i(d)  =  p0  =  fej  ■ 

•  P4  :  {dx(j,  k)}  Short  Range  Dependent  if  n  >  H  + 1/2. 

\2qk  —  2J  k'\  — +  +00,  |Cov dx(j,  k)dx(j',  k')\  <  D\2qk  —  2J  k'\2<yH~N\ 

N  >  1  and  Dilation 


=*>  Idealisation  :  dx(j,k)  Independent  Variables  . 

=*  Interpretations:  X(t)  =  Ekax(J,k)‘pj,k(t)+Ej=i,...,j,kdxti, *0^00- 

=>  Implications:  E\dx(j,k)\q  =  E\dx{0,k)\q2jqH  \fq/E\dx(0,k)\q  <  oc. 
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Wavelets  and  Long  Range  Dependence 


•  Spectral  Analysis  : 

Let  X  be  a  2nd  Order  stationary  process, 

Let  VI;  be  the  FT  of  ip  with  central  frequency  u0  and  bandwith  Au0. 

E|dx(i,  k)\2  =  frx(v)\*(Vv)\dv 

~  2_jTx(2_jVo)  within  bandwith  2~jAiy0. 


Let  X  be  Long  Range  Dependent 


-  Power  Law:  Tx(v)  =  cf \v 

-  Power  Law:  E\dx(j,k)\2  - 


—  a 


0  <  a  <  1, 


V 


->  0 


C2^°l~1\  j  -►  +oo, 


•  {dx(j,  k)}  Short  Range  Dependent  if  X  >  a  - l. 

|2Jfc  —  2J  k'\  —y  -poo,  Cov dx (j.  k)dx(j\  k')\  <  D\2^k  —  2J  k' 

N  >  1  and  Dilation 


cl  —  1  —  27V 
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2nd  Order  wavelet  Statistical  Analysis 


Abry,  Gongalves,  Flandrin 


Principles: 

-  Ideas  :  PI  =>•  E| dx(j,  k) |2  =  C22j2H 

=>  log2  E\dx{j,k)\2  =  j2H  +  j3q, 

-  Problems:  Estimate  IE | dx(j,  k)\2  from  a  Single  Finite  Length  Observation  ? 

-  Solution  :  P2  et  P3  =>  Statistical  Averages  =>  Time  Averages, 

S2U)  =  (l/n:)Z:U\dxti,k)\2 

Log-Scale  Diagrams:  iog2 s2(j)  vs iog2 t  =  j 
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2nd  Order  Wavelet-based  Statistical 
Analysis  for  Self  -Similarity 
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2nd  Order  Wavelet-based  Statistical 
Analysis  for  Long  Range  Dependence 


Octave  j 
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Wavelets  and  2nd-order  Scaling:  Estimation 


Dyadic  Grid  (Discrete  Wavelet  Transform)!  cij  =  2 =  k23 


SCALE 


TIME 


Structure  Function  (Time  Average): 

Yj  =  (51°g25'2(2j)  =  |log2(l /nj)J2kU  \dx(j,k) |2 


Definition 


Y.  versus  log2  2j  =  j, 


H  = 


3=3 1  3  3 


WHERE  E,  jwj  =  1,  E,  Wj  =  0,  WITH  Wj  = 


AND  p  =  0,1,2,  Bp  =  EjiP/ai!  aj  ARBITRARY  NUMBERS. 


What  Are  the  Performance  of  Such  an  Estimator  ? 

When  applied  to  a  Self-Similar,  or  LRD  Process 
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Wavelets  and  2nd-order  Scaling:  Estimation 


Abry,  Gongalves,  Flandrin,  Abry,  Veitch 

•  ASSUME:  -  i)X  Gaussian, 

-  ii)  Idealisation:  exact  independence. 

•  Bias  :  Elog2  S2(j)  =  log2  ES2(j)  +  T'(nj/2)  -  \og2(nj/2\. 

'' - - v - y 

9j 

=>  EH  =  H  +  \  wjgj. 

•  Variance:  -Var  H  = 

—  min  Var  H  =^>  aj  oc  Var  log2  S2U) 

—  Var  log2S,2(j)  —  C/rij  ~  2^C /n, 

=>■  VAR  H  ~  ((log2e)2(]T\M;22J))  /n, 

=>  Analytical  (approximate)  Confidence  Interval 

(DOES  NOT  DEPEND  ON  UNKNOWN  H). 

•  Actual  Performances  :  negligible  biais,  extremely  close  to  mle. 

•  Conceptual  and  Practical  Simplicity  :  matlab  code  Available. 
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Wav.  and  2nd-order  Scaling:  Robustness 


Superimposed  Trends 

Y(t)  =  X(t)  +  T(t)  =>  dY{j,k)  =  dx{j,k )  +  dT(j,k) 


-  If  T(t)  Polynomial  of  degree  P,  then  dT  =  0  when  N  >  P, 

-  If  T(t)  smooth  trend,  then  the  dT  decrease  as  N  increases. 


30 


25 


20 


10 


-101 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - L- 1 

2000  4000  6000  8000  10000  12000  14000  16000 


Vary  N  ! 
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Wav.  and  2nd-order  Scaling:  Robustness 


Superimposed  Trends  -  Ethernet  Data  (Veitch,  Abry) 


x  105 
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Wav.  and  2nd-order  Scaling:  Robustness 


Constancy  along  time  of  Scaling  laws  (Veitch,  Abry) 


pAug 


pAug 


Time 

(j1=7,j2=12) 
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Self-Similarity 


•  Self-Similarity: 

E\dx(j,k)\q  =  Cq(2*yH 

-  Power  Laws, 

-  V2J  (for  all  scales), 

-  Vq/E\dx(j,  k)\q  <  cx), 

-  A  single  parameter  H 

-  Additive  Structure. 


? 


[33] 


Beyond  Self-Similarity 


•  Self-Similarity: 

E|  dx(j,k)\*  =  Cq(2?yH 

-  Power  Laws, 

-  V2J  (for  all  scales), 

-  Wq/E\dx(j,  k) \q  <  oo, 

-  A  single  parameter  H 

-  Additive  Structure. 

•  Multi  Fractal 

E\dxU,k)\i  =  Cq(V)^ 

-  Power  Laws, 

-  V2J  <  L,  (for  fine  scales  only,  in  the  limit  2j  — >  0,) 

-  Vg? 

-  A  whole  collection  of  scaling  parameter  C(g) 

-  Multiplicative  Structure. 


? 
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Modelling  Tool  2:  Multiplicative  Cascades 


•  Definition: 

-  Split  Dyadic  Intervals  lj:k  into  two, 

-  I.I.D.  Multipliers  Wjik 

~  Qj(t)  =  U{(j,k):l<j<J,teIjjk}^j,k, 


Implications: 

-  Local  Holder  Exponent, 

-  MultiFractal  Sample  Paths,  MultiFractal  Spectrum  D(h ) 


-  Cascades,  Multiplicative  Structure, 

-  Efc  [}/a tik+aT° x(u)duY  =  <EMCg, 


Fine  Scales  a  — >  0, 


-  Multiple  Exponents  Cg, 

-  No  Characteristic  Scale, 

-  Cg  =  -  log2  EWq,  Non  Linear  in  q. 
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Modelling  Tool  2:  Multiplicative  Cascades 


Yaglom,  Mandelbrot 

Mandelbrot’s 
Cascade  (CMC) 

-  1 1 D  W, 

-  Dyadic  Grid, 


Barral,  Mandelbrot  Schmmitt  et  al., 


Compound  Poisson 
Cascade  (CPC) 

-  1 1 D  W, 

-  Point  Process, 


Bacry  et  al.,  Chainais  et  al. 
Infinitely  Divisible 
Cascade  (IDC) 

-  Continuous  Infinitely 

-  Divisible  Measure  M, 


Qr(t)  =  n  Wjfkj 
<p(q)  =  -  log2 


n  Wj^n 

=  -q(  1  -  IE W)  +  1  -  IE Wq, 
A(t)  =  limr^0  Jo  Qr(u)du, 


exp  f  dM(t' ,  r) , 

=  p(q)  ~  qp(  1), 


For  a  Range  of  qs,  IE |  A(t  +  ar0)  -  A(t) \q  = 

Resolution  Depth  <  Scale  <  Integral  Scale,  a7 


jm 


5 

r  <  a  <  aM  =  L. 


Cq  |  a 
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MultiFractal  Processes 


Density:  Qr(t)  =  n  wJk 

E  (i  ftt+aT°  Qr{u)du )  =  cqa^\ 

Measure:  A(t)  =  limr^0  fo  Qr(u)du, 

IE| A(t  +  ar0)  -  A(t)\q  =  cq\a\q+v{q\ 


Fractional  Brownian  Motion 
in  Multifractal  Time: 

VH(t)  =  BH(A(t)), 

E| VH(t  +  ar0)  -  VH(t)\q  =  cq\a\qH+^qH\ 

MF  Process  t 

1.6 
1.4 

MultiFractal  Random  Walk:  : 

Y„(t)  =  f*  Qr(s)dB„(s),  Z 

E| YB(t  +  ar0)  -  Y„(t) |«  =  c,|a| - 


Matlab  Synthesis  Routines  :  Chainais,  Abry 
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Higher-Order  Wavelet  Statistical  Analysis 


Principles  : 

-  Ideas  :  PI  =>■  El dx(j,  k) \q  =  E\dx(0,  k)\q2K<i 

=>•  log2  IE | dx(j,  k)\q  =  jCq  +  (3q, 

-  Problems:  Estimate  IE | dx(j,  k) \q  from  a  Single  Finite  Length  Observation  ? 

-  Solution  :  P2  et  P3  =>  Statistical  Averages  =>  Time  Averages, 

SqU)  =  ^/nj)YZ1\dxU,k)\q 

Log-Scale  Diagrams!  \og2 sq(j)  vs iog22j  =  j 
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LogScale  Diagrams  -  Multi  Fractal  Proc. 


Octave  Octave  Octave  Octave 
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Wav.  and  higher-order  Scaling:  Estimation 


•  Dyadic  Grid  (discrete  wavelet  transform):  a,  =  2j,tjtk  =  k2\ 

SCALE 

•  •  • 

•  •  •  •  • 


TIME 

•  Structure  Functions  (Time  Average): 

sqti)  =  0-/nj)YZi\dxU,k)\* 

•  Definition: 

Yj,q,n  =  log2  Sn( 2\  q-  fo)  VERSUS  log2  =  j, 

C(<hn)  =  Sj=21  wj,qYj,q,n  ■ 

Non  Weigthed:  a3  =  cste 

•  What  Are  the  Performance  of  Such  Estimators  ? 

When  applied  to  MultiFractal  Processes 
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'EST  FOR  THE  FlNITENESS  OF  MOMENTS 


Gonqalves,  Riedi 


Theorem  : 


Let  X  BE  A  RV  WITH  CHARACTERISTIC  FUNCTION  x(s)  :  =  IE  exp{isX}. 

If  H^x  :=  sup{a  >  0  :  |9ftx(s)  —  Pa(s)  \  <  C\s\a}, 

IS  THE  LOCAL  HOLDER  REGULARITY  OF  AT  THE  ORIGIN,  THEN 

IE  | X | g  <  +c X)  Vg  <  <?+  AND  H&x  <  q+  <  \H^X\  +  1. 


Estimator  : 

{Xfc}fc=i....,„,  n  I.I.D  RVS,  SET 

n 

W (a)  :=  n~1  ^^(a.Xfc) 

k= 1 


WITH  ^  A  REAL  AND  SEMI-DEFINITE 

Fourier  transform  of  a 

SUFFICIENTLY  REGULAR  WAVELET  ^ . 

Then 


a,  log  |  W(o)| 

rLmx  =  lim  sup - 

a — >0+  & 
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Estimating  the  Partition  Function  Support 
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Methodology 


•  Numerical  Synthesis  of  Processes: 

-  Accumulate  nbreal  numerical  replications  with  length  n  samples. 

•  Apply  Scaling  Exponents  Estimators: 

-  Compute  C(g,  n)^  for  each  replication, 

-  Average  over  repl.  to  obtain  the  statistical  performance  of  C(g,  n) 

•  Asymptotic  Behaviours: 

-  The  cascade  depth  increases  for  a  given  number  of  Integral  Scales. 
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Methodology 


•  Numerical  Synthesis  of  Processes: 

-  Accumulate  nbreal  numerical  replications  with  length  n  samples. 

•  Apply  Scaling  Exponents  Estimators: 

-  Compute  £(<?,  n)^  for  each  replication, 

-  Average  over  repl.  to  obtain  the  statistical  performance  of  C(g,  n ) 

•  Asymptotic  Behaviours: 

-  The  cascade  depth  increases  for  a  given  number  of  Integral  Scales. 

-  The  number  of  Integral  Scales  increases  for  a  given  cascade  depth, 
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Linearisation  Effect:  C(g) 

Lashermes,  Abry,  Chainais 

CPC  Qr  EI(l)  CPC  VH  Eli  1(3) 


q  >  Qo ,  n)  =  a0  +  /30q,  qQ,  ot0,  Po  ARE  RV. 
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Linearisation  Effect:  Legendre  Transform 


D(h)  =  d  +  MIN q(qh  -  £(<?)),  (d  EUCLIDIEN  DIMENSION  OF  SPACE). 


CPC  Qr  EI{1) 


CPC  VH  Eli  1(3) 


Accumulation  points  :  D0(h0 ),  with  Dq  =  d  -  aQ,  hQ  =  /30, 

D0 ,  hQ  ARE  RV. 
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Lin.  Effect:  Asymptotic  behaviours 


Given  resolution,  Increasing  Number  of  Integral  Scales, 


qo 


9  10  11  12  13  14  15  16  17 

log2(n) 


•  Given  Number  of  Integral  Scales,  Increasing  Resolution,  q0 


13  14  15  16  10  11  12  13  14  15  16  10  11  12  13 

log2(n)  log2(n)  log2(n) 


[48] 


Linearisation  Effect:  Conjecture 


•  Critical  Points: 


(  Df  =  0, 

J  D(hf)  =  0, 

[  hf  =  (, dC(q)/dq)q=q± . 


Results: 


El 


C(q,  n)  =  d  -  D~  +  h~q 
C(q,  n) 

C(q,  n)  =  d-  D+  +  h+q 


EIIk.HI  : 


C(<2,  n) 

C(«>  n)  =  d  -  D+  +  h+q 


d-  D~  +  h~q,  q  <  q~ , 

C(«),  q*  <  q  <  q t> 

d  -  D+  +  h+q,  q+  <  q. 

->  C(q),  0  <  <3  <  q+, 

d  -  D+  +  h+q,  q+  <  q. 
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Linearisation  Effect:  Comments 


When  does  the  Linearisation  Effect  EXIST  ? 

-  FOR  ALL  TYPES  OF  CASCADES:  CMC,  CPC,  IDC, 

-  FOR  ALL  TYPES  OF  PROCESSES:  Qr,  A,  VH,  YH, 

-  FOR  ALL  NUMBERS  OF  VANISHING  MOMENTS:  N  >  1, 

-  FOR  ALL  MRA-BASED  ESTIMATORS:  WAVELETS,  INCREMENTS,  AGGREGATION, 

-  CAN  BE  WORKED  OUT  FOR  q  <  0, 

-  EXTENDS  TO  DIMENSION  HIGHER  THAN  d  >  1. 
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Extension:  Standard  WT  versus  WTMM  (1/3). 


Standard  vs  WTMM 
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Extension:  2D  Multiplicative  Cascade  (2/3). 
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Extension:  3D  Multiplicative  Cascade  (3/3). 


3D  CMC  (log  Normal),  El(1)  compared  to  a  2D  slice. 
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Linearisation  Effect:  Comments 


When  does  the  Linearisation  Effect  EXIST  ? 

-  FOR  ALL  TYPES  OF  CASCADES:  CMC,  CPC,  IDC, 

-  FOR  ALL  TYPES  OF  PROCESSES:  Qr,  A,  VH,  YH, 

-  FOR  ALL  NUMBERS  OF  VANISHING  MOMENTS:  N  >  1, 

-  FOR  ALL  MRA-BASED  ESTIMATORS:  WAVELETS,  INCREMENTS,  AGGREGATION, 

-  CAN  BE  WORKED  OUT  FOR  q  <  0, 

-  EXTENDS  TO  DIMENSION  HIGHER  THAN  d  >  1. 

What  the  Linearisation  Effect  IS  NOT: 

-  A  LOW  PERFORMANCE  ESTIMATION  EFFECT. 

-  A  FINITE  SIZE  EFFECT  :  THE  CRITICAL  PARAMETERS  DO  NOT  DEPEND  ON  n, 

BE  IT  THE  NUMBER  OF  INTEGRAL  SCALES, 

OR  THE  DEPTH  (OR  RESOLUTION)  OF  THE  CASCADES. 

-  A  FINITENESS  OF  MOMENTS  EFFECT, 

-  q~  <  0  <  1  <  <?+,  q  -  1  +  ip(q)  =  0, 

-  Qc  <  q*  <  0  <  1  <  <it  <qt> 

What  the  Linearisation  Effect  Might  Be: 

-  Multiplicative  Martingales  ? 

-  Ossiander,  Waymire  00,  Kahane,  Peyriere  75,  Barral,  Mandelbrot  02. 
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Linearisation  Effect:  Picture 


Two  POWER-LAWS,  TWO  FUNCTIONS  OF  q\ 

Bare  Cascade:  EQr{t)q  =  rv(q),  q  €  n. 


Dressed  cascade: 

ETQo(t,  a;  /30)q 
E TQo(t,  a;  f30)q 


cM(q\ 


=  oo, 


Q  €  [qc  ,  <?+], 

ELSE, 


with: 


C(«) 

C(g) 

C(g) 


*  5 


1  +  qh 

<p(q), 

1  +  qht 


q  G 
g  G 


[g; 

[g; 


*  5 


fc  ’  g* 

Q  €  [g+,  g+ 


]. 

], 

]• 


•  Confusion  between  <^(g)  and  C(g): 

-  Multiplicative  Cascade:  ip(q),q  e  TZ, 

-  Scaling  Exponents:  <(g),  q  €  [gj,  g+]. 
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Linearisation  Effect:  Sketched  Views 


Moments 


oo 


i  i 

I  <  OO  I 

-I - 1- 


C  T 

-q — 


m 


<  OO  I  oo 


I 


E  A  (t) 

T 


q  _ 


qc  q*  -101  q*  q* 


Estimated  C(q>^qh 


m 


1+qh 


+ 


Estimated  £(q,_n) 
Ell  &  EMI 


«q) 


1+qh 


+ 


- 1 
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Linearisation  Effect:  Impacts  and  Importance 


Consequences:  Recast  the  Usual  Goals  : 

-  Estimate  the  Integral  Scale  and  the  Resolution  of  the  Cascade, 
=>  i.e. ,  Find  a  Scaling  Range  [am,  aM\ 

-  Estimate  the  Critical  Parameters  Df,  hf,  qf, 

-  Estimate  the  CO?)  for  q  e  [q~ ,  g+], 

— >  Visit  B.  Lashermes’s  Poster. 


Importance  of  the  Linearisation  Effect: 

-  Discrimination  of  MF  Models  based  on  C(g,  n), 

-  Discrimination  between  monoFractal  and  MultiFractal, 
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Negative  Values  of  qs 


Difficulties  ? 

_  Ti  ' 

-  Finiteness  ?  Sq(j)  =  (l /nj)  Y.kU  I dx(j,  k)\q  <  oo? 

-  Numerical  Instability  ?  dx(j ,  fc)  ~  0  — >  fc)lg  =  oo 

-  Theory  ?  Full  MultiFractal  Spectrum 


q  q 


Solutions  ? 
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Negative  Values  of  qs  -  Solution  1 


Aggregation!  Tx(a,t)  =  f*+aT°  x(u)du 


Applies  only  to  positive  data  (Measure) 
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Iog0(a)  logja) 


Negative  Values  of  qs  -  Solution  2 


WT  Modulus  Maxima  (arneodo  etal.) 

Lx(a,tk )  =  SUPa/<a|Tx(a/,  tk(a'))\ 


WTMM 

N\ 


0  100  200  300  400  500  600  700  800  900  1000 

Skeleton  of  the  Wavelet  Transform 


0  100  200  300  400  500  600  700  800  900  1000 

Time 


Computationally  Expensive 
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Negative  Values  of  qs  -  Solution  3 


Wavelet  Leaders:  (jaffard  etal.) 


o  o 


o  o 


dx(j ,  k)  — >  Lx{j ,  fe)  =  2  J  ) 


o  o 


o  o 


o  o 


Computationally  Efficient  and  Excellent  Statistical  Performance 


[61] 


Beyond  Power  Laws 


•  Self-Similarity: 

E\dx(j,  k)\q  =  Cq(  2j)qH  =  Cq  exp(qH  In  2j) 

-  Power  Laws, 

-  V2J  (FOR  ALL  SCALES), 

-  Vq/E\dx(j,  k)\q  <  oo, 

-  A  SINGLE  PARAMETER  H 

-  Additive  Structure. 

•  MultiFractal 

E\dx(j,k)\q  =  Cq{ 2*)C(,)  =  Cqex p(C(q)  In  2J) 

-  Power  Laws, 

-  V2j  <  L,  (FOR  FINE  SCALES  ONLY,  IN  THE  LIMIT  2j  — >  0,) 

-  Vg? 

-  A  WHOLE  COLLECTION  OF  SCALING  PARAMETER  £(g) 

-  Multiplicative  Structure. 


•  Beyond  Power  Laws  :  Warped  Inf.  Div.  Cascades 

E| dxU,  k)\q  =  Cq(yyH  =  Cq  exp(qH  In  2-’) 
E\dx(j,k)\q  =  Cq{ 2*)c(,)  =  Cq  exp(£(q)  In  V) 
E|dx0',fc)|9=  =  C,exp(C(q)n(2J)) 

— >  Visit  Pierre  Chainais’s  Poster 
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Conclusions  and  References 

Analysing  Scaling  in  Data  ? 

-  Think  wavelet 

—  EFFICIENCY, 

—  PRACTICAL  AND  CONCEPTUAL  ADEQUATION  AND  SIMPLICITY, 
—  ROBUSTNESS  AGAINST  NON  STATIONARITIES, 

—  EASY  TO  USE,  LOW  COAST,  REAL  TIME  ON  LINE. 

Modelling  Scaling  in  Data  ? 

-  Think  Self  Similarity  versus  Multiplicative  Cascades, 

-  AND  POSSIBLY  ADD  LONG  MEMORY. 

-  also  Scaling  may  not  be  power  laws 


References  and  Resources,  Visit  : 

—  perso.ens-lyon.fr/patrice.abry 

—  inrialpes.fr/is2/~pgoncalv 

—  www.cubinlab.ee.mu.oz.au/~darryl 

—  f  raclab 

—  www.isima.fr/~chainais 
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